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Abstract 

The properties of the two-dimensional exactly solvable Lieb and 
Baxter models in the critical region are considered based on the ther- 
modynamic method of investigation of a one-component system criti- 
cal state. From the point of view of the thermodynamic stability the 
behaviour of adiabatic and isodynamic parameters for these models 
is analyzed and the types of their critical behaviour are determined. 
The reasons for the violation of the scaling law hypothesis and the 
universality hypothesis for the models are clarified. 

Key words: scaling law hypotheses; universality hypotheses; sta- 
bility coefficients; critical state 

1 Introduction 

The description of the behaviour of thermodynamic parameters near the crit- 
ical points is one of the basic problems of the critical state theory. Direct 
statistical calculations connected with the evaluation of the partition func- 
tion of real systems are unavailable at present because of the impossibility 
of accounting exactly for the interactions and, moreover, for the fluctuations 
which are large near the critical point. So, solving the problem by the meth- 
ods of statistical physics one considers either the simplest models, for which 
the partition function can be evaluated exactly, or an approximate solution 
of the problem. 
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At the first approach the exactly solvable two-dimensional models (the 
Ising, Lieb, Baxter models and others [T] forming the most valuable posses- 
sion of statistical mechanics) are of great importance. The second approach 
is connected mainly with the examination of the asymptotic behaviour of 
thermodynamic parameters near the critical points, as well as with the de- 
velopment of the scaling law hypothesis, the universality hypothesis and the 
renormalization group approximation in various variants and has appreciably 
succeeded. Indeed, the large class of real systems and models satisfies the 
scaling law and the universality hypotheses. The existence of real systems 
and exactly solvable two-dimensional models, for which these hypotheses are 
violated is also remarkable. The six-vertex ferroelectric Lieb model and the 
eight-vertex Baxter model [T] are such examples. 

Our aim is the examination of the critical properties of these models 
based on the thermodynamic method of investigation of the critical state 
[2]-[l] which is developed on the first principles without any hypotheses. 

The method is based on the constructive critical state definition and the 
critical state stability conditions. The method describes a variety of critical 
state nature manifestations. The violation of the scaling law and universality 
hypotheses in the Lieb and Baxter models is explained just by this variety. 

The Lieb and Baxter models give a reasonable fit to real ferroelectrics 
(antiferroelectrics) and ferromagnets (antiferromagnets). So, the application 
of the thermodynamic method to them could be interesting for the 

development of the critical state theory. 

2 The thermodynamic method of investiga- 
tion of the critical state 

Let us consider the basic theses of the thermodynamic method and the ter- 
minology. The critical state definition, which considers both the properties 
of homogeneous and heterogeneous system can be written in the form [2]^I1]: 

Here X is the generalized thermodynamic force, x is the conjugated ther- 
modynamic variable (the external parameter of a system), K^. is the critical 



2 



slope of a phase equilibrium curve. Eq. ([T]) has non-trivial solutions, if the 
condition 
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(2) 



is fulfilled all over the spinodal. It coincides with the well-known critical 
state condition D = 0, where D is the stability determinant of the system 
[H |6] . According to the terminology of Refs. O [6] the parameters concerned 

, , , , . . , , fdT\ fdT\ fdX\ 

under the constant thermodynamic variables, I I , I I , I I , 

are the adiabatic parameters (AP's); the parameters concerned under the 



constant thermodynamic forces 
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namic parameters (IP's) 
the adiabatic stability coefficients (ASC's); whereas 



dx ) 

The parameters 7— 



'dX\ , . , 

— — , are the isody- 

are called the isodynamic stability coefficients (ISC'c). The stability coeffi- 
cients are related to the fluctuations of the external parameters of the system 
(the first and the second Gibbs lemmas) which infinitely increase near the 
critical point. 

The definition ([T]) describes the critical state by means of the AP's. The 
solution of the homogeneous linear equations ([1]) is the critical slope K,. 
which distinguishes the critical point on the spinodal. It is the fundamental 
characteristic of the critical state and it can be expressed via the ASC's: 



dx 



sign 



dx) s. 




(3) 



This definition, being combined with the critical state stability conditions, 
leads to the existence of four alternative types of the critical behaviour of 
thermodynamic systems [2]^[1]- The behaviour type is defined by the value 
of one ASC and Kc- 

The behaviour of the whole set of the stability characteristics of the sys- 
tem (the AP's and IP's) is determined for each type. The fourth type of the 
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critical behaviour is the most interesting and the most "fluctuating" one. In 
this case it is necessary to consider the differential equations of higher orders. 
Then the solution is realized by several possibilities |2]-|1]. The case of two 
or even three phase equilibrium curves converging at the critical point is of 
special interest. Such a point has not yet been found experimentally, but 
in this paper we demonstrate that the critical point of the ferroelectric Lieb 
model has just this feature. 

3 The ferroelectric 6-vertex Lieb model 

There are a lot of crystals with the hydrogen bonds in the nature [T]. The 
ions in such crystals (with the coordination number four) must obey the ice 
rule. The bonds between atoms via hydrogen ions form the electric dipoles. 
So, it is convenient to represent them as the arrows on the bond curves. 
These arrows are directed to that end of the bond which is occupied by the 
ion. There are only six such configurations of arrows, therefore the ice models 
are sometimes called the six-vertex models. The partition function of such a 
system is defined by the expression 



where the summation should be carried out over all the configurations of 
the hydrogen ions allowed by the ice rule, is the energy of i-tjpe vertex 
configuration and Ui is the number of i-type vertices in the lattice. 

There are three sorts of the ice models which have been solved by E. H. 
Lieb [71 |S]. One of them, considered in this paper, can describe KH2PO4 
(KDP), the crystal with hydrogen bonds, which is characterized by the coor- 
dination number four and orders ferroelectrically at low temperatures under 
the appropriate choice of ei, £2, • • • , ^e- For the square lattice this choice is 



In the ground state all the arrows are directed either up and to the right or 
down and to the left. Both these states are typical for the ordered ferroelec- 
tric. 

The expression for the free energy per lattice point in the presence of the 
nonzero external field is given by 



Z = exp[-(ni£i + n2e2 + . . . + neeQ)/kT], 



(4) 



€1=62 = 0, £3 = £4 = £5 = £6 > 0. 



(5) 
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Figure 1: The phase diagram of Lieb model [T] 



where P is the electric polarization |T], and A = — 0.2122064/cTc; k is Boltz- 
mann constant. The critical equation of state is expressed in the form 



E 



P 



k{T-%) 
sign(£^) 



if \E\ < k{T - Te) 
otherwise. 



(7) 



It corresponds to the phase diagram in Fig. [T] 

It is necessary to emphasize that the ice model allows the investigation 

on the basis of the thermodynamic method. In this case the temperature T 

and the electric intensity E stand for the generalized thermodynamic forces. 

The conjugated generalized thermodynamic variables are the entropy S and 

the electric polarization P. Thus, the adiabatic parameters for the given 

, fdE\ , , . , 
and I -r^ I , and the isodynamic parameters 



model are 



dT 
dS 



[dP 



dP, 

As T ^ 



the free energy per lattice 



dT\ fdT\ , fdE 

point coincides with expression ([6]), and as T — )■ the free energy equals 
simply to Si — EP. Consequently, the heat capacity is finite in the subcritical 
region and the critical exponent is a' = 0. Both the phases are quite ordered 
and then differ from each other only by a direction of the electric polarization 
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vector (P = ±1). This corresponds to the second critical behaviour type 
according to the thermodynamic classification of critical behaviour types of 

"no ) ~ 7^ 7^ \ "no ) ~ ^' Thus, 

the critical slope of the equilibrium curve of the phases / and // equals to 
zero, Kc = 0. 

As we can see from Eq. ([6]), in the supercritical region (T — t- T+) the 
heat capacity diverges as I — — — - ) , i.e. the thermic ASC is I — — 



T, ' \dS 



Cy — — — . Let us approach to the critical point from the supercritical 

region along the curve of the first-kind phase transition / — III and // — ///. 
It is known that at least one of the jumps AP, AS" must exist along these 
curves. I.e., on the transition curve 

AP = P,-P,„ = l-^^^0. (8) 

At the critical point AP = 0. 

The entropy jump can be determined from the known behaviour of the 
heat capacity. For the phase / we have a' = 0, i.e. Cp = const. Conse- 
quently, the entropy of the phase I is Sj = CilnT + const. For the phase 
/// we have a = 1/2, i.e. Sm = C2^/Tc{T — Tc) + const. Then, for the 
jump, we have 



AS = Si- Siii = Ci In T-C2 VTc{T - + const. (9) 

At the critical point AS* = const 7^ {0, 00}. Such a behaviour of the entropy 
is connected with the divergence of the heat capacity in the supercritical 
region. 

The analogous results can be obtained for phases // — /// as well. 
For the equilibrium line I — II we have AP = 2, AS* = 0. At the critical 
point AP = 0. 

Thus, the found values of the jumps correspond to the results of papers 
[2]-[l], and the point is critical for the phase equilibrium line / — //, for 
the line / — /// and for the line II — III. So, the point C in the phase 
diagram (Fig. [1]) is the point of the convergence of three phase equilibrium 
lines. The possibility of such a point has been predicted in papers [2]-|l]. 
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Let us analyze the behaviour of the whole set of the system stability 
characteristics (the AP's and the IP's). The relations between the adiabatic 
and isodynamic parameters exist: 



dS L\dp L~ [dS L\dP ~~\dP AdP 



(10) 



Using Eqs. ([7]) and ( JTOjl . we can obtain the following expressions for the AP's 
and the IP's: 
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/dT\ T-Tc 
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where C = -§ = 6.2831908 ■ f . The critical slope equals kP = kP for 

the line I - III and kP = -kP for the line // - ///. 

As it follows from Eq. ( fTTl) . at T — t- all the thermodynamic stability 
characteristics tend to zero: 

dT\ fdE\ fdT\ 

^ ^0, Uttt ^0, ^0, 



dS) p ' "' \dP) s ' \dP) s ' 
dT\ ^ (dE\ ^ (dT\ 

Os)e^'^ 

According to the critical behaviour classification [3] at -f^c 7^ {0, oo} and 
ASC's— )■ we have the fourth type of the critical behaviour, and two phase 
equilibrium lines with different critical slopes Kc = ±kP converge at the 
critical point. This behaviour type is the most fluctuating one (the fluctu- 
ations of the energy and polarization {AHy, (AP)^ — )■ oo). Approaching 
to the critical point from the subcritical region (along the phase equilibrium 
line I — II with the slope Kc = 0) the second type of critical behaviour is 
realized (the fluctuations of the energy {AHy is finite and the fluctuations 
of the polarization (AP)^ — )■ oo). 

As it is known, stability characteristics are inversely proportional to fluc- 
tuations of external parameters of the system. At the continuous transitions 
[6j D and the SCs pass finite minima, that corresponds to the growth of 
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fluctuations. The locus of these minima is curve of supercritical transitions 
(the lowered stability curve or quasispinodal). The limit case of these contin- 
uous transitions, when fluctuations in the system are at the high and D and 
the SC's pass zero minima, is the critical state. The critical point is also the 
limit point of some first-kind transition (the limit point of phase equilibrium 
curve). If the phase equilibrium curve and curve of supercritical transitions 
pass into each other continuously, i.e. the slopes of these curves are the same, 
then the tricritical point is observed, where three phases become identical: 
two subvritical phases and supercritical one. 

On the quasispinodal the next condition is fulfilled [9]: 



Using results (fTTl) to find the determinant of stability for Lieb model, 
and investigating where condition (fT2l) is fulfilled, we obtain = 0. This is 
equation of quasispinodal for ferroelectric Lieb model. The resulting phase 
diagram is shown in Fig. [2j So, the maximal growth of fiuctuations is 
observed under zero electric field. The critical slope of the subcritical phase 
equilibrium curve is ii"c = 0. It means, that for this model it is realized the 
case of continuous passage of the equilibrium curve into the lowered stability 
curve because of the same critical slopes. 

Thus, the violation of the scaling law hypothesis in the Lieb model can 
be explained by the fact that the model corresponds to two different critical 
behaviour types: at T — )• T+ the second type and at T~ the fourth type is 
fulfilled. Besides, the critical point of the Lieb model is the critical point of a 
special type with the convergence of three phase equilibrium lines. Moreover, 
the equilibrium curve continuously passes into the lowered stability curve. 

4 8-vertex Baxter model 

The eight-vertex Baxter model is a generalization of the six-vertex Lieb 
model. The ice models as the models of the critical phenomena have some un- 
usual properties: the ferroelectric state at these models is frozen (i.e. there is 




(12) 



or, equivalently. 




8 




Figure 2: The quasispinodal for Lieb model (the gray hne) 

complete ordering even at the non-zero temperature); the critical behaviour 
of antiferroelectrics is characterized by the more complicated law instead of 
a simple power dependence of (T — Tc). 

The first of these unusual properties is connected with the ice structure. 
In the case of an unlimited lattice with the ferroelectric ordering the infinite 
energy is needed for a deformation. So, the deformation gives an infinitesimal 
contribution to the partition function [Tj. 

The following generalization of the ice-type models was proposed [TO]^ 

• there is only one arrow on each square lattice edge; 

• the configurations with an even number of arrows getting in (and get- 
ting off) each vertex are allowed only; 

• eight possible configurations of the arrows to a vertex exist. 

The formation of j-type vertex needs the energy ej (where j = 1, 8). 

For such a model the partition function is given by (jl]) where the sum- 
mation is performed over the eight vertex configurations. 

Thus, besides the first six vertices coinciding with the Lieb model there 
are another two new vertices for which all the arrows get either in a vertex 
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or off a vertex. Now the finite energy is needed for the local deformation of 
the lattice state (e.g. for reversing of all the arrows which are lying on the 
square side), in which all the arrows are directed up or to the right. So the 
ferroelectric state is ordered not completely. 

As it was mentioned above, the Baxter model is fitted to describe the crit- 
ical phenomena in ferroelectrics (antiferroelectrics). The eight- vertex model 
can be considered also as two Ising models with the nearest neighbours in- 
teraction (each model is on its sublattice). These sublattices are connected 
by means of the four-spin interplay. In this case the model corresponds to 
ferromagnets. 

The Baxter model has the exact solution only in the absence of an external 
field. The critical exponents of this model equal [T] 



71 

a = a = 2 



„ 7r-/i 7r + /i 37r + /i 

4/i 2/i IT — jj, 

Here the index e denotes the electric exponents. The exponents (3, 7 and 5 
are related to ferromagnet. The exponent a is the same both for the ferro- 
magnet and for the ferroelectric, /i is the interaction parameter, it takes a 
value from (0, vr). Thus, as we can see, the critical exponents depend on the 
interaction parameter continuously. This fact is in the contrary to the uni- 
versality hypothesis. This result distinguishes the Baxter model among the 
others two-dimensional exactly solvable models. Taking this into the consid- 
eration, one should expect that the type of the critical behaviour according 
to the thermodynamic classification and the value of the critical slope change 
depending on the interaction parameter. Let us show this. 

4.1 The ferromagnetic Baxter model 

In the case of ferromagnet the adiabatic stability coefficients get the following 
asymptotic form: 

TT TtT 

dT\ 2-- /dH\ ^ 
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where t 



T-Tr 



It is necessary to note that in the absence of the ex- 



ternal field the behaviour of the isodynamic parameters coincides with the 

vr 

behaviour of the adiabatic parameters. When < /i ^ — the exponent a is 
negative, the exponent 7 takes a positive value, i.e. 



dH 



dT 
dM 



s 



Kr, = 



TT 1 StT 

and the second type of critical behaviour is fulfilled. At — < u < 

the fourth type of critical behaviour is realized, the exponent a increases 

14 7 14 

(0 < a < — ) and the exponent 7 decreases (- > 7 > — ), where a < 7. 




Figure 3: The temperature dependence for reduced thermic coefficient of 
stability 
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Figure 4: The temperature dependence for reduced magnetic coefficient of 
stability 



(dH \ ( dH \ 
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tend to zero faster than other parameters. The value of the critical slope is 

ISvr 

Kc = 0. The case fi = corresponds also to the fourth type of critical 

16 
14 

behaviour, but a = 7 = — and all the parameters tend to zero according 

15 

ISvr 

to the same law, the critical slope is Kc 7^ {0, 00}. At < < vr the 

14 14 
fourth type of critical behaviour is also fulfilled, — < a < 1 and — > 7 > 

15 15 

T ( dT 

- and everywhere a > 7. All the parameters tend to zero, but I —— 

and I 7-— ) tend to zero faster than other parameters. The value of the 



dS 



H 



critical slope is K^. = 00. The corresponding plots of ASC's for various /i are 
presented in Figs. [3], HI 

TT 

Thus, the performed analysis has determined that at < ^ — the crit- 
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ical behaviour of the Baxter model corresponds to the second type according 
to the thermodynamic classification [2]^[1] with = 0, and at — < /i < vr 

it corresponds to the fourth type which is fulfilled by three possibilities for 
the critical slope {K^. = 0, Kcj^ {0, oo}, Kc = oo) depending on the value of 
/i varying within the mentioned interval. 



4.2 The ferroelectric Baxter model 

In the case of the ferroelectric Baxter model the stability coefficients can be 
written in the form: 

71 vr + /i 

dT\ 2-- /g^ 



dS J p ' V-^^ 



s 



TT 

At < ^ — , as in the previous case, a is negative and 7 is positive. So 
dT \ ( dE \ ( &T \ 

— ) 7^ 0, ( — ) = ^ ( — ) = 0, i^c = and the second type of 



dS ) p' ' \dP ) s \dP J s 

7T 

critical behaviour is fulfilled. At — < < vr the exponent a takes positive 

3 

values < a < 1, but a is less than 7, - > 7 > 1 and the fourth type of 

critical behaviour with J^c = is realized. 

It is interesting to emphasize the fact that for real ferromagnets and 
ferroelectrics the critical behaviour types are also the second and the fourth 
ones. 



5 Conclusion 

Thus, in the paper the consideration of the thermodynamic stability of the 
Lieb and Baxter models by the method of Ref. [2]-[l] has been performed. 
The asymptotic expressions for the whole set of the stability characteristics 
are determined. The reasons for the violation of the scaling law and univer- 
sality hypotheses in the given models are clarified. So, we determine that the 
second and the fourth type of critical behaviour is fulfilled in the subcriti- 
cal and in the supercritical region of the Lieb model, correspondingly. The 
violation of the scaling law hypothesis in the ferroelectric Lieb model can 
be explained just by the difference of the behaviour types. It has been also 
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ascertained that three phase equihbrium hnes with different critical slopes 
converge at the critical point of the model. A possibility of the existence of 
such a type of the critical point has been predicted in papers |2]-|1]. The 
equation of quasispinodal is obtained and it is shown that the equilibrium 
curve continuously passes into the lowered stability curve in this model. 

In the Baxter model the fulfillment of the second and the fourth type of 
critical behaviour also occurs, moreover, the fourth type is represented by 
three possibilities — with three different critical slopes of the phase equilib- 
rium line. The reason for the violation of the universality hypothesis is that 
each of the mentioned types (the second type, the fourth type with Kc = 0, 
the fourth type with 7^ {0, 00} and the fourth type with = 00) is 
connected either to the certain value or the continuous range of the inter- 
action parameter /i. It is interesting to emphasize that in each model while 
one hypothesis is violated the another one is nevertheless holds. In addition, 
the special case of the eight-vertex Baxter model, in which the universality 
hypothesis is violated, is the Lieb model [fi = 0), where the universality hy- 
pothesis is satisfied, but the scaling law hypothesis is violated, and the Ising 
model (/i = — ), where both hypotheses are fulfilled. Therefore, the abilities 

of the thermodynamic method of investigation of the one-component system 
critical state have been illustrated by the example of the above-mentioned 
models and the global reasons for the violation of the scaling law and uni- 
versality hypotheses concerned with the variety of the critical state nature 
manifestation are revealed. 
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